Spin-1/2 Heisenberg antiferromagnet on an anisotropic kagome lattice 
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We use the coupled cluster method to study the zero-temperature properties of an extended 
two-dimensional Heisenberg antiferromagnet formed from spin-1/2 moments on an infinite spatially 
anisotropic kagome lattice of corner-sharing isosceles triangles, with nearest-neighbor bonds only. 
The bonds have exchange constants Ji > along two of the three lattice directions and J2 = «Ji > 
along the third. In the classical limit the ground-state (GS) phase for k < 1/2 has collinear 
fcrrimagnetic (Neel') order where the J2-coupled chain spins are ferromagnetically ordered in one 
direction with the remaining spins aligned in the opposite direction, while for k > 1/2 there exists 
an infinite GS family of canted ferrimagnetic spin states, which are energetically degenerate. For 
the spin-1/2 case we find that quantum analogs of both these classical states continue to exist as 
stable GS phases in some regions of the anisotropy parameter ft, namely for < ft < tt ci for the 
Neel' state and for (at least part of) the region ft > ft C2 for the canted phase. However, they 
are now separated by a paramagnetic phase without either sort of magnetic order in the region 
ft ci < ft < ft C2 , which includes the isotropic kagome point k = 1 where the stable GS phase is now 
believed to be a topological (Z2) spin liquid. Our best numerical estimates are ft ci = 0.515 ± 0.015 
and ft C2 = 1.82 ± 0.03. 

PACS numbers: 75.10.Jm, 75.30.Gw, 75.40.-s, 75.50.Ee 



I. INTRODUCTION 

Low-dimensional quantum magnets, especially those 
defined on regular two-dimensional (2D) lattices, have 
been the subject of intense study in recent years (and 
see, e.g., Refs. [l|Hfor recent reviews). In particular it is 
known that highly frustrated 2D quantum antifcrromag- 
nets display a bewildcringly rich panoply of ground-state 
(GS) phases, which often have no classical counterparts. 
Examples include various valence-bond crystalline and 
spin-liquid phases. Among the parameters that deter- 
mine the zero-temperature (T = 0) phase diagram of 
such systems are the dimensionality and structure (e.g., 
the coordination number) of the crystallographic lattice 
on whose sites the magnetic ions are situated, the spin 
quantum number s of the ions, and the type and range of 
the magnetic bonds between the ions that often compete 
for differing forms of order, thereby leading to frustra- 
tion. We have learned too that the quantum versions of 
classical models that have massively degenerate ground 
states, especially those with a nonzero (T = 0) GS en- 
tropy, are prime candidates for systems with novel GS 
phases. 

Among all such candidate spin-lattice systems, there- 



fore, those with periodic arrays of vertex-sharing struc- 
tures, each of which is itself magnetically frustrated, oc- 
cupy a special niche. These include the three-dimensional 
(3D) pyrochlore lattice of vertex-sharing tetrahedra and 
the 2D kagome lattice of corner-sharing triangles. Of 
these, the spin-1/2 Heisenberg antiferromagnet (HAF) 
on the 2D kagome lattice has been the subject of inten- 
sive study in recent years Even after several decades 
of research the nature of the GS phase of the spin-1/2 
HAF on the spatially isotropic kagome lattice has re- 
mained uncertain until very recently. Various outcomes, 
ranging from states with magnetic order to valence-bond 
solids or quantum spin liquids of different types, have 
been proposed. 

Although many such studies agree on the finding that 
the GS phase lacks magnetic long-range order (LRO), 
there has remained uncertainty over its precise character. 
For example, some studies have favored a gapless criti- 
cal spin liquid of one type or another, while others have 
favored one or other valence-bond crystals with an abun- 
dance of low-lying spin-singlet excited states. Only in the 
last year or so has compelling numerical evidence been 
provided^ due to advances in the density-matrix renor- 
malization group (DMRG) technique, that the ground 
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state is both gapped and is without any signal of ei- 
ther valence-bond or magnetic order at the largest finite- 
size systems that could be studied. In the past few 
months further convincing evidence has come from dif- 
ferent large-scale DMRG studie s 49 ' 50 that this GS phase 
is indeed a topological (Z2) spin liquid, as we discuss 
further in Sec. fVl below when we discuss our own results. 

Theoretical interest in spin- 1/2 kagome HAFs height- 
ened considerably in the last few years with the dis- 
covery of several candidate materials for their experi- 
mental realizations. Chronologically, the first promis- 
ing such candidate was the mineral herbertsmithite (also 
known as Zn-paratacamite), 7-Cu3Zn(OH)6Cl2r^ r — for 
which it has been shown that the spin- 1/2 Cu 2+ ions 
are antiferromagnetically coupled and lie on the vertices 
of well separated and structurally undistorted kagome- 
lattice planes. Although the underlying kagome planes 
in herbertsmithite appear to be essentially structurally 
perfect, there does appear to be an appreciable amount 
of antisite disorder due to a mixing of the spin- 1/2 Cu 2+ 
ions and the diamagnetic Zn 2+ ions between the Cu and 
Zn sites. This disorder acts to introduce a coupling be- 
tween the kagome planes, thereby effectively destroying 
the local 2D nature of the system. Thus, while herbert- 
smithite is structurally perfect, these impurities, together 
with a spin-orbit coupling that can be modelled by a 
Dzyaloshinskii-Moriya interaction with a non-negligible 
strength parameter, act to complicate the comparison 
of theory with experiment. As a consequence herbert- 
smithite has lost some of its initial promise as an almost 
perfect, spin- 1/2, isotropic kagome HAF. 

A more recently discovered candidate for that role is 
another member of the atacamite family, namely the 
polymorph kapellasite, a-Cu3Zn(OH) 6 Ci2r^— of her- 
bertsmithite. Although they share the same chemical 
composition, the two minerals have a different crystal- 
lographic structure. Interestingly, however, they both 
display distinct kagome structures, although in differ- 
ent ways. Thus, in kapellasite the spin-1/2 kagome lat- 
tice is obtained by the regular doping of a 2D triangu- 
lar Cu 2+ metal-site sublattice with diamagnetic Zn 2+ 
ions. By contrast, in herbertsmithite it is obtained by 
a similar diamagnetic dilution with Zn 2+ ions of the 3D 
pyrochlorc-likc sublattice. It is asserted 5 -! for kapellasite 
that while the Cu/Zn mixing leads to some intralayer dis- 
order within the kagome planes, it cannot induce any ap- 
preciable intcrlayer coupling, unlike in herbertsmithite. 

It should be noted, however, that a theoretical elec- 
tronic study using density functional theory (DFT) 
within the local density approximation^ 5 - of both 
the material kapellasite and its relative haydeeite, 
Cu3Mg(OH)6Cl2, has revealed significant non-NN ex- 
change coupling strengths, especially those correspond- 
ing to bonds across the diagonals of the hexagons on the 
kagome lattice. Furthermore, recent high-temperature 
series expansion fits to the measured DC magnetic sus- 
ceptibility, XDc(r), as a function of temperature T, for 
kapellasite^ seem to give a nearest-neighbor (NN) ex- 



change interaction ( Ji) on the kagome planes that is fer- 
romagnetic in nature (i.e., J\ < 0), with the overall anti- 
ferromagnetic behavior of the material explained by large 
positive further-neighbor interactions. 

A spatially anisotropic version of the spin-1/2 
kagome HAF has also been suggested to have 
been realized experimentally in the minerals vor- 
borthite, Cu 3 V20 7 (OH)2-2H 2 0^^ and vesignieite, 
BaCu 3 (V0 4 )2(OH)2.£ 2 -£^ In both of these materials 
the Cu sites form a slightly distorted kagome network, re- 
sulting in two inequivalent Cu sites per triangle, namely 
one Cul site and two Cu2 sites per triangle. For ex- 
ample, volborthite has a monoclinic distortion that de- 
forms the equilateral triangles of the isotropic kagome 
network into isosceles triangles. In this material the dif- 
ference in the Cul-Cu2 and Cu2-Cu2 bond lengths is 
about 3%. In turn, theoretical modelling of the thermo- 
dynamic properties then leads to a suggested difference 
between two of the NN (Cul-Cu2) exchange constants 
(Ji) and the third (Cu2-Cu2) one (J[) on each trian- 
gle of the kagome lattice. In volborthite this magnetic 
anisotropy is around 20%. The anisotropy is much less 
pronounced in vesignieite where the difference in bond 
length is less than 0.1%, and the material is closer to 
being structurally isotropic. Another recently discovered 
spin-1/2 deformed kagome-lattice antiferromagnet is the 
material Rb 2 Cu 3 SnFi 2 ^-^ 

Although both volborthite and vesignieite have a re- 
duced symmetry compared with the structurally perfect 
herbertsmithite, they do offer some advantages. As noted 
above, this latter compound shows antisite disorder with 
up to 10% of the magnetic Cu 2+ ions exchanged by Zn 2+ 
ions, thereby leading to a weak interlayer magnetic cou- 
pling between the kagome planes as well as magnetic 
vacancies within them. By contrast, in both volbor- 
thite and vesignieite their intermediate layers between 
the kagome planes contain V 5+ ions, and hence antisite 
disorder of the Cu ions is prevented. Although vesig- 
nieite is much less anisotropic than volborthite it suffers 
in practice, like herbertsmithite, from low sample qual- 
ity. One of the main experimental advantages of studying 
the more anisotropic volborthite over cither of herbert- 
smithite or vesignieite is that it is much easier to prepare 
with fewer impurities. Nevertheless, we note that the na- 
ture of the magnetic couplings in this material has been 
questioned in a recent study^ where it is pointed out 
that the local environments of the two inequivalent types 
of Cu sites differ essentially in important ways. DFT is 
then used to show that volborthite should not be mod- 
eled as an anisotropic J\-J[ kagome-lattice HAF, but 
rather as a J^—J^—Ji model, more reminiscent of coupled 
frustrated chains, in which two-thirds of the kagome sites 
(viz., the Cu2 sites) are considered as J[-J2 chains (i.e., 
with ferromagnetic NN exchange, J[ < and frustra- 
tion induced by next-nearest-neighbor (NNN) exchange, 
J2 > 0), and with the chains coupled via NN exchange 
bonds of strength J\ between the Cu2 and remaining Cul 
sites. 



3 




FIG. 1: (Color online) (a) The interpolating kagome-square model with (black) solid bonds — Ji, (green) dashed bonds - - J[, 
and (red) dash-dot bonds - ■ - J2, showing the canted state; and the equivalent anisotropic kagome model when J[ — 0, showing 
(b) the coplanar ferrimagnetic canted state and (c) the collinear ferrimagnetic semi-striped state. In all cases the (blue) arrows 
represent spins located on lattice sites •. 



Although it is thus often very uncertain as to whether 
a given real material does or does not provide an experi- 
mental realization of a particular theoretical model, such 
as the kagomc-latticc HAFs considered here, there is still 
much to be gained by a systematic theoretical comparison 
between such models. It is of particular interest in such 
comparisons to use, wherever possible, the same theoret- 
ical technique. Among the relatively few widely applica- 
ble and systematically improvable (within a well-defined 
hierarchical approximation scheme) such tools is the cou- 
pled cluster method (CCM).— — Our intention here is to 
use the CCM to further the study of the HAF, with NN 
interactions only, on an anisotropic kagome lattice) 7 - 7 - - — 

By now the CCM has been used to study a huge 
number of quantum s pin- lattice problems (see, e.g., 
Refs . fl^l^l4ai7347M8l-ll03 and references cited therein). 
Among these, of particular interest here are applica- 
tions of the CCM to the frustrated spin-1/2 Ji~J 2 HAF 
model on the square lattic o 85 ' 92 ~ 94 i 98 with NN bonds (of 
strength Ji > 0) competing with NNN bonds (of strength 
Ji > 0), and various models related to it by removal 
of some of the NNN J 2 bonds. When half of the J 2 
bonds are removed these include an interpolating square- 
triangle HAF (or spatially anisotropic triangular HAF)^ 
the Union Jack lattice model^ and the anisotropic pla- 
nar pyrochlore (or checkerboard) HAF (alao known as 
the crossed chain model) 

A further modification of the original square-lattice 
J1-J2 model is now to remove another half of the J 2 
bonds in such a way as to leave half of the fundamen- 
tal square plaquettes with one Ji bond and the other 
half with none. One way of doing this in a regular 
fashion results in the Shastry-Suthcrland modeli^ in 
which no Ji bonds meet at any lattice site and every 
site is five-connected (by four NN J\ bonds and one J2 
bond). The CCM has also been successfully applied to 
this modeh 89 ' 97 

Another similar such model, of particular interest here, 
arises from removing alternate diagonal lines of J2 bonds 
from the interpolating square-triangle HAF (which itself 



arises from the square-lattice J\~Ji HAF model by re- 
moving all of the diagonal lines of Ji bonds in the same 
direction). It differs from the Shastry-Sutherland model 
primarily in that the square lattice now breaks into two 
square-sublattices of A sites and B sites, respectively, 
such that the A sites are all six-connected (by four NN 
J\ bonds and two NNN Ji bonds), while the B sites are 
all four-connected (by four NN J\ bonds only). A par- 
ticularly relevant generalization of the model for present 
purposes arises from introducing an additional anisotropy 
in the NN bonds such that along alternating rows and 
columns the NN bonds are allowed to have the strength 
J{, as shown in Fig. [If a). Clearly, when J[ = the 
model reduces to the anisotropic kagome-lattice HAF, 
shown equivalently in Fig. HJb), which is the subject of 
the present paper. 



The spin-1/2 interpolating kagome-square model de- 
scribed above and shown in Fig.[TJa) was studied by us in 
an earlier paper, using the CCMj2& In that paper we were 
mainly interested in the quantum phase transition line in 
the J[-Ji plane (with J\ = 1) in the model between the 
two quasiclassical states with antiferromagnetic Neel or- 
der and ferrimagnetic canted order. By contrast, in the 
present paper we focus purely on the J[ = case corre- 
sponding to the anisotropic kagome-lattice HAF that has 
been suggested as one possibility to describe the mag- 
netic properties of volborthite and vesignieite, as dis- 
cussed above. Our main aim is not so much to shed light 
on the structure of the paramagnetic GS phase of the 
isotropic spin-1/2 kagome-lattice HAF, but to determine 
the boundaries of this phase as the anisotropy is varied. 



After first describing the model in Sec. [Hi we apply 
the CCM to investigate its GS properties. The CCM 
is itself described briefly in Sec. IIII1 and our results are 
then presented in Sec. HVl We conclude in Sec. [V] with a 
discussion of the results. 



4 



II. THE MODEL 

As mentioned in Sec. HI in a previous paper— we consid- 
ered a depleted (and anisotropic) variant of the archety- 
pal and much-studied J\-J 2 model in which three- 
quarters of the J2 bonds are removed from it in the pat- 
tern shown in Fig. [TJa). The square-lattice representa- 
tion of the model shown in Fig. HJa) contains the two 
square sublattices of A sites and B sites respectively, and 
each of these in turn contains the two square sublattices 
of Ai and A2 sites, and Bi and B2 sites respectively, as 
shown. It is very illuminating to consider the anisotropic 
variant (viz., the interpolating kagome-square model or 
Ji-J[-J2 model) in which half of the J\ bonds are al- 
lowed to have the strengths J[ > along alternating rows 
and columns, as shown in Fig.[lja). All of the bonds join- 
ing sites i and j are of standard Heisenberg type, i.e., pro- 
portional to Sj-Sj, where the operators = [s%,s v kl s%) 
are the quantum spin operators on lattice site k, with 
s| = s(s + 1) and s = | for the quantum case considered 
here. 

The spin-1/2 HAF's on the 2D kagome and square lat- 
tices are represented respectively by the limiting cases 
{.h = J 2 ,J{ = 0} and {Ji = J[, J 2 = 0}. The limit- 
ing case when { J\ — J[ — 0; J 2 > 0} represents a set of 
uncoupled one-dimensional (ID) HAF chains. The case 
J[ = with J 2 7^ Ji represents a spatially anisotropic 
kagome HAF considered recently by other authors^ - — 
especially in the quasi-lD limit where J2/J1 li 80 ' 82 It 
is this latter model where J{ = that is considered here 
(although, for technical reasons, we note that we actually 
set J[ to be a small positive value, henceforth chosen to 
be J[ = 10~ 5 ). Our model is thus equivalently shown in 
the kagome-lattice geometry of Fig. [ljb). It thus com- 
prises parallel chains of spins on A-sitcs coupled in a NN 
fashion along the chains by bonds of strength J2, with 
the chains then cross-linked via B2 sites with NN bonds 
of strength Ji . The totality of sites arranged in a regu- 
lar kagome lattice on which each component triangle is 
comprised of two J\ bonds and one J 2 bond. 

The Hamiltonian of the resulting anisotropic kagome- 
lattice model is thus 

H = Ji ^2 s i ■ s j + <h Si ' S J ' W 

(hi) 

iGAj'eB 2 ieAiJeA 2 

where the sum on runs over all NN pairs (of the sort 
specified in each sum), counting each bond once and once 
only. Henceforth we consider the model where the spins 
on all lattice sites have spin quantum number s = ^ , and 
where both types of bonds are antiferromagnetic in na- 
ture (i.e., J 1 > 0, J 2 > 0) and hence act to frustrate one 
another. With no loss of generality we may then choose 
the energy scale by setting J\ = 1. We are interested in 
the infinite-lattice limit, Nk = — > 00, where Nk is 
the number of sites on the kagome lattice and N is the 
number of sites on the square lattice of Fig. [TJa) before 
the (jN) Bi sites have been removed. 



Considered as a classical model (corresponding to the 
case where the spin quantum number s — > 00) the in- 
terpolating kagome-square model of Fig. (Ha) (i.e., with 
J[ 7^ 0) has only two GS phases separated by a con- 
tinuous (second-order) phase transition at J2 = Jf 1 = 
^(Ji + J[). For J 2 < J 2 the system is Neel-ordered 
on the square lattice, while for J 2 > J 2 the system 
has noncollincar (but coplanar) canted order as shown 
in Fig.[]Ta), in which the spins on each of the Ai and the 
A2 sites are canted respectively at angles (it =f <j>) with 
respect to those on the B sublattice, all of the latter of 
which point in the same direction. The lowest-energy 
state in the canted phase is obtained with <fi = (f> c \ = 
cos -1 ^ 1 /^). The Neel state, for J 2 < Jf 1 , simply cor- 
responds to the case <f> c \ = 0. 

For the case of the anisotropic kagome-lattice HAF 
considered here (i.e., with J[ = 0), the classical (s — > 00) 
ground states are those spin configurations that satisfy 
the condition that for each elementary triangular plaque- 
tte of the kagome lattice in Fig. [1Tb) the energy is mini- 
mized. If we take the angle <fi to be such that the middle 
spin (viz., that on a B2 site) of a given triangular plaque- 
ttc forms angles (ir±(f>) with the other two (chain) spins of 
the same plaquette (viz., those on A2 and Ai sites respec- 
tively), the total energy of the lattice, for classical spins 
of length s, is E = |A^ A -s 2 [2Ji cos(tt - 0) + J 2 cos(2</>)]. 

For J 2 < \J\ this energy is minimized with = 0, 
and the classical GS phase is thus collinear and unique, 
with the spins (on the A sites) along the J2-bond chains 
aligned in one direction and the remaining spins (on 
the B2 sites) on the kagome lattice aligned in the op- 
posite direction. As a convenient shorthand notation 
we henceforth refer to this collinear fcrrimagnetic state 
as the Neel' state. Indeed, the Neel' state of the 
anisotropic kagome-lattice HAF is precisely equivalent to 
the Neel state of the interpolating kagome-square model 
of Fig. QTa) (i.e., before the removal of the Bi sites in 
the limiting case J[ = 0) from which it is derived. The 
total spin of this classical collinear Neel' ferrimagnetic 
state is thus Stot = ^Nks where each spin has magni- 
tude s. In terms of the saturation magnetization (i.e., in 
the ferromagnetic state with all spins aligned in the same 
direction), M sat = Nks, the total magnetization in this 
collinear ferrimagnetic state is M tot = -|M sat . For the 
quantum case the Marshall-Lieb-Mattis theore m 106 ! 107 
may also be used to show, for the limiting case J 2 = 
only, that the exact ground state has the same value 
Stot = ^Nks of the total spin as its classical counter- 
part. 

By contrast, for J 2 > the classical GS en- 

ergy is minimized with the canting angle <f) = <p c \ = 
cos_1 (^)- We expect that coplanar canted states will 
then be favored by either thermal or quantum fluctua- 
tions, and henceforth we only consider coplanar states 
from among this degenerate manifold that includes non- 
coplanar states. The total magnetization of this canted 
ferrimagnetic state is M tot = | (2 cost/) — 1)Nks = 
|(x _ 1 )^sat- The ground state of the HAF on the 
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isotropic kagome lattice (i.e., with J2 = Ji) falls in this 
regime, and has the canting angle <p = § demanded by 
symmetry. Only for this case docs the total classical 
magnetization vanish, Af tot = 0. The classical ensem- 
ble of degenerate coplanar states is now characterized by 
two variables for each triangular plaquette, namely the 
angle <f>, and the two- valued chirality variable x = ±1 
that defines the direction (anticlockwise or clockwise) in 
which the spins turn as one transverses the plaquette in 
the positive (anticlockwise) direction. For a given value 
of J2 > \ J 1 the different degenerate canted states arise 
from the various possible ways to assign positive or nega- 
tive chiralities to the triangular plaquettes of the lattice. 
(Appendix A of Ref. [ll gives a good description of the 
constraints that these chiralities need to satisfy.) 

We note that the well-known q = and -\/3 x -y/3 
states of the (isotropic) kagome-lattice HAF are the spe- 
cial cases, respectively, where all of the chiralities are the 
same, and where basic triangular plaquettes joined by a 
vertex have opposite chiralities. The q = state is one in 
which the spins on each of the three sublattices (of Ai, 
A2, and B2 sites respectively) are parallel to one another, 
and make an angle of 120° with the spins on the other 
two sublattices, while the x state contains nine 
sublattices and is obtained by deleting | of the sites (viz., 
the Bi sites) of the ordered spins of a triangular lattice to 
form the kagome lattice (and see also, e.g., Refs. H043| 
for further details). Clearly, the limiting case J[ —> of 
the classical ground state of the interpolating kagome- 
square model shown in Fig. QJa) is just the state shown 
in Fig.[IJb), which corresponds to the q = state for the 
isotropic (J2 = J\) kagome-lattice HAF. 

The HAF on the isotropic kagome lattice (i.e., with 
J2 = J\) is especially interesting since for this case, with 
<f> = |, the number O of degenerate classical spin con- 
figurations grows exponentially with the number Nk of 
spins, so that even at zero temperature the system has a 
nonzero value of the entropy per spin. A previous high- 
order CCM study^ of the isotropic kagome-lattice HAF 
showed that for the extreme quantum case, s = h, the 
q = state is energetically favored over the \/3 x \/3 
state, while for any s > i the \/3 x \/3 state is selected 
over the q = state. For both the \/3 x \/3 and the 
q = states it was further found that the magnetic or- 
der is strongly suppressed by quantum fluctuations. In 
particular, the order parameter (viz., the avearage local 
on-site magnetization or sublattice magnetization) Mk 
was found to vanish for both s — i and s = 1, while 
nonzero values for Mk were found for s = |, 2, |, and 
3. 

By contrast, for the anisotropic case (with J2 7^ Ji), 
the classical degeneracy Q has been shownZ^ to grow ex- 
ponentially with \/Nk [i-e., £1 oc exp^^/Nx)}, so that 
the GS entropy per spin vanishes in the thermodynamic 
limit. Clearly, since in the limit J 2 — > J\ the anisotropic 
model approaches the isotropic model, the anisotropic 
model must have an appropriately large number of low- 



lying excited states that become degenerate with the 
ground state in the isotropic limit, J 2 — > J\- 

We also note that as J2 — > 00 the classical canting an- 
gle (j) c \ — > J 71 ") an d the spins on the A sublattice chains be- 
come antifcrromagnctically ordered, as is expected, and 
these spins are orientated at 90° to those on the B2 sub- 
lattice. For the particular ordering in Fig. [Tfb) , which 
arises from that in Fig. [Ha) in the limit J[ — > 0, the 
spins on the B2 sublattice are themselves parallel and 
hence fcrromagnctically ordered. Of course there is com- 
plete degeneracy at the classical level in this decouplcd- 
chain limit (i.e., when J 2 — > 00) between all states for 
which the relative ordering directions for spins on the A 
and B2 sublattices are arbitrary. In the same limit the 
quantum spin- 1/2 problem considered here should also 
comprise decoupled antiferromagnetic chains on the A- 
sublattice sites. We expect that this degeneracy in rela- 
tive orientation might be lifted by quantum fluctuations 
by the well-known phenomenon of order-by- disorder.-}^ 
Since it is also true that quantum fluctuations generally 
favor collinear ordering, a preferred state is thus likely to 
be the so-called ferrimagnetic semi-striped state shown 
in Fig. QTc) where the A sublattice is now Neel-ordered 
in the same direction as the B2 sublattice is fcrromagnct- 
ically ordered. (Note that in the square-lattice geometry 
of Fig. [Ua) and where the Bi sites and the J[ bonds 
are retained, alternate rows (and columns) are thus fer- 
romagnetically and antiferromagnetically ordered in the 
same direction in the semi-striped state, which is the ori- 
gin of its name.) 



III. COUPLED CLUSTER METHOD 

We now apply the CCM (see, e.g., Refs. frMlE and ref- 
erences cited therein) to the spin-1/2 anisotropic kagome- 
lattice HAF discussed in Sec. |H] above. At a very gen- 
eral level the method provides one of the most versa- 
tile techniques now available in quantum many-body the- 
ory. At attainable levels of computational implementa- 
tion it has been shown to provide some of the most ac- 
curate results ever obtained for a large number of quan- 
tum many-body systems in quantum chemistry, as well 
as in condensed matter, atomic, molecular, nuclear, and 
subnuclear physics i 72 ' 73 More specifically for present pur- 
poses, it has been very successfully applied by now to 
a large number of systems of interest in quan tum mag- 
netism (see, e.g., Refs. HHHlillTl-tzMil-llol and refer- 
ences cited therein), as we have already noted in Sec. |TJ 

The method of applying the CCM to quantum mag- 
nets has been described in detail elsewhere (see, e.g., 
Refs. [HlzlziHii] and references cited therein). It 
relies on building multispin correlations on top of a suit- 
ably chosen, normalized, GS model (or reference) state 
|$) in a systematic hierarchy of approximations that we 
described below. The reference state |$) is required only 
to be a fiducial vector for the system in the sense that all 
possible states of the system can be described in terms 
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of it as a linear combination of states obtained from it 
by acting on it with members of some suitably chosen 
complete set of mutually commuting multispin creation 
operators, Cj = (Cj - )*. In this way |$) acts as a gen- 
eralized vacuum state with respect to the set of opera- 
tors {C+}. It is often chosen as a classical ground state 
of the model under investigation, and for the present 
anisotropic kagomc-latticc HAF we use mainly the canted 
(coplanar) fcrrimagnctic state shown in Fig. [ljb) (includ- 
ing the collinear Neel' state which is its limiting form 
when the canting angle <f> — > 0), although we also discuss 
briefly in Sec. [V] the use of the semi-striped (collinear) 
ferrimagnetic state shown in Fig. [TJc) as a CCM model 
state. 

Once the set {|$),C/} has been suitably chosen, the 
exact GS ket and bra wave functions of the system are 
parametrized within the CCM, in terms of them, in the 
exponential forms that are the hallmark of the method, 

|*}=c s |$), S = J2SiCl, (2) 
(*| = ($|£c- s , 5=1 + 53&C7, (3 ) 

where we define C + =1. It is clear from Eqs. @ 
and (|3|) that the normalization has been chosen so that 
(*|*) = ($|$) = ($|$) = 1. The complete set of GS 
CCM correlation coefficients {Si, Si} (VI ^ 0) is then 
obtained by the requirement that the states (\E'| and |$) 
obey the GS Schrodinger equations, (^\H = E(^\ and 
H\^) = E\^>), respectively. The resulting equations, 

($|C7e~ s iJc s |$) =0, VI ^0, (4) 
($|5e- s [iJ,C+]c s |$) =0, V7V0. (5) 

may, completely cquivalently, be derived from the re- 
quirement that the GS energy functional 

H=(^\H\iSf), (6) 

be stationary with respect to variations in all members 
of the set {Si,S~i;I ^ 0}. 

Once the CCM correlation coefficients have been found 
from solving Eqs. (jj) and §Q it is easy to see that the 
GS energy E is given purely in terms of the ket-state co- 
efficients {Si} as E = (<S>\e~ s He s \&). Clearly, however, 
for a more general operator O, the evaluation of its GS 
expectation value, O = (^f\0\^f), requires knowledge of 
the set of bra-state correlation coefficients {Si} as well 
as of the corresponding set of ket-state coefficients {Si}. 

It is very convenient in practice to perform a rotation 
of the local spin axes of each of the spins in the system 
(i.e., we define a different set of spin axes on every lattice 
site) such that all spins in the reference state align along 
the negative z axes of the local coordinates. In practice 
this simply means that the Hamiltonian has to be rc- 
cxprcssed in terms of these local axes for each choice of 



reference state used. The big advantage of so doing is 
that in these local coordinates we have 

|$> = UW--->, Cf = 8+8+ ■■■8+ , n=l,2,3,..., 

(7) 

where = + is%, the indices k n denote arbitrary 
lattice sites, and the components of the spin operators 
are defined in the local rotated coordinate frames. We 
note that for spins of quantum number s, each site in- 
dex k n in each multispin configuration set-index / = 
{fci, k2, ■ ■ ■ k n } in Eq. (|7|) can be repeated up to a maxi- 
mum of 2s times. Thus, for the present case, s = |, all 
individual spin-site indices in each set-index / are differ- 
ent. 

The magnetic order parameter (viz., the average 
local on-site magnetization) is now given by M = 

~W Y^,iLi(^\ s t W' wnere s i is expressed in the local spin 
coordinates defined above, and N(— > 00) is the number 
of lattice sites. We denote by M here the order param- 
eter defined for the general interpolating kagome-square 
model of Fig. QJa) (i.e., before the jN number of Bi 
sites have been removed to yield the anisotropic kagome 
model). The corresponding order parameter, Mk , for the 
anisotropic kagome-lattice HAF considered here is simi- 
larly given by M K = - Y,?=i (*l s l I*) > where the sum 
is taken over only the kagome-lattice sites. Thus, in the 
limiting case J[ = considered here of the anisotropic 
kagome-square model, the spins on the non-kagome Bi 
sites are frozen to have their spins aligned exactly along 
their local negative z axis. Hence, for the anisotropic 
kagome-lattice limit (i.e., when J[ = 0) of the interpolat- 
ing kagome-square model we have the simple relation 



The parametrizations of Eqs. @ and yield, in prin- 
ciple, the exact GS eigenstate when the complete sets 
of multispin creation and destruction operators, C+ and 
Cj , respectively, is retained. In practice, of course, it 
is necessary to make approximations by truncating the 
complete set of multispin configuration set-indices {/}. 
In that case the results for physical quantities such as 
the GS energy E and order parameter M will naturally 
depend both on the particular truncation (i.e., on the 
configurations specified by the set-indices I that are re- 
tained), as well as on the specific choice of model state 
to which those multispin configurations are referred. We 
note, however, that the CCM always exactly obeys the 
Goldstonc linked-cluster theorem at every such level of 
approximation^ Thus, the CCM approach always yields 
results directly in the thermodynamic limit, N — > 00, 
from the outset, and no finite-size scaling is required. 
The results at all levels of approximation are guaranteed 
to be size-extensive. 

We note, however, that although the CCM is fully (bi- 
)variational as discussed above in connection with the 
stationarity of H in Eq. ©, it does not lead to strict 
upper bounds for the energy due to the lack of explicit 
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hcrmiticity in the paramctrizations of the GS kct and bra 
wave functions in Eqs. ^ and ([3]). This minor drawback 
is more than compensated for by the fact that the CCM 
parametrizations exactly obey the important Hellmann- 
Feynman theorem at all levels of approximation^ 

For the present spin- 1/2 model we employ the so-called 
LSUBm (or lattice-animal-bascd subsystem) approxima- 
tion scheme to truncate the expansions of S and S in 
Eqs. ([2]) and ([3]). In this very widely tested scheme one 
includes from the full set of multispin configurations spec- 
ified by the set-indices / in Eqs. ©, ([3]), and ([7]) only 
those involving m or fewer correlated spins in all arrange- 
ments (or lattice animals in the language of graph the- 
ory) which span a range of no more than m contiguous 
lattice sites. In this context a set of sites is defined to 
be contiguous if every site has at least one other in the 
set as a nearest neighbor, and where it is clearly neces- 
sary to include a definition in the geometry (or better, 
topology) of the lattice of which pairs of sites are con- 
sidered to be NN pairs. For example, we choose here to 
work in the triangular-lattice geometry of the anisotropic 
kagome-square model in which the B sublattice sites of 
Fig. (Ha) are defined to have four NN sites joined to them 
by either Ji bonds or J[ bonds, and the A sublattice sites 
are defined to have the six NN sites joined to them by 
Ji, J[, or J 2 bonds. If we had chosen instead to work 
in the square-lattice geometry every site would have four 
NN sites. The former triangular-lattice geometry leads 
in the limiting case when J[ = to the natural kagomc- 
lattice geometry of Fig. [T^b) in which any two of the Ai, 
A2, and B2 sites forming the basic triangular plaquettes 
are considered to be NN pairs. Clearly, this would not be 
the case in the latter square-lattice geometry. For more 
details of the LSUBm scheme the reader is referred, for 
example, to Refs. [tHe! 

Clearly, the LSUBm truncations scheme provides a 
fully systematic approximation hierarchy in the sense 
that each time the truncation index m is increased more 
of the Hilbert space is sampled, such that asm->oo the 
approximation becomes exact. Although, as we have al- 
ready indicated, we never need to perform any finite-size 
scaling, since all CCM approximations are automatically 
performed from the outset in the infinite-lattice limit, 
Nk —> 00, where Nk is the number of lattice sites, we 
do still need as a last step in a CCM calculation to ex- 
trapolate to the exact m — > 00 limit in the LSUBm trun- 
cation index m, at which the complete (infinite) Hilbert 
space is reached. By now there is a great deal of experi- 
ence available regarding how one should extrapolate the 
GS energy per site eo(m) = E(m)/NK and the magnetic 
order parameter M(m). 

Thus, for the GS energy per spin, eo(m), we use the 
well-tested empiri cal scal ing ansatz (and see, e.g., Refs. 
[86ll89ll90ll92l49 J96ll9§ 4l^olflol . 

eo(m) = «o + aim~ 2 + ci2m~ 4 . (9) 

For the GS magnetic order parameter different extrap- 
olation rules have been used depending on the features 



of the spin-lattice system at hand. For highly frustrated 
spin-lattice systems, a well-tested rule (and see, e.g., 
Refs. [92ll9ll98l fl0Cl - [l0l is 

M(m) = b Q + b im ~ 1/2 + b 2 m- 3/2 . (10) 

An alternative rule that is useful, for example, for situ- 
ations when there is some frustration present, but when 
it is not too large, i o 89 ' 104 

M(m) = c + c x m^ v , (11) 

which, as an advantage, leaves the leading exponent open 
for determination. The disadvantage of this scheme is 
clearly that it involves only the leading two terms in the 
asymptotic expansion, compared to the (inherently more 
accurate) three terms used in Eq. (| 10[) . Clearly, however, 
when the exponent v in Eq. (jll[) is found by fitting to the 
CCM LSUBm results to be close to the value 0.5, as is 
usually found for highly frustrated systems, we can then 
revert to the more accurate form of Eq. (JTUJ) . 

In the present paper we present extrapolated results 
based on LSUBm data sets with m = {2,4,6,8}. To 
check the robustness of the extrapolation rules and to es- 
timate the associated error bars in the extrapolations we 
have also performed extrapolations using LSUBm data 
sets with m = {2,4,6} and m = {4,6,8}. We find in 
general that the extrapolated results from all three sets 
of data are very similar. 

The number of independent fundamental LSUBm clus- 
ters (i.e., those that are inequivalent under the symme- 
tries of the Hamiltonian and of the model state) that 
are retained in the expansions of Eqs. © and ([3]) for 
the CCM correlation operators S and S increases rapidly 
with the truncation index m. For example, the number 
of such fundamental clusters for the canted model state 
of the interpolating kagome-square model of Fig. Ufa) 
is 201481 at the LSUB8 level of approximation in the 
triangular-lattice geometry used here where J2 bonds are 
considered to join NN pairs, and this is the highest level 
for the present model that we have been able to attain 
with available computing power. In order to solve the 
corresponding coupled sets of CCM bra- and ket-state 
equations we use an efficient parallelized CCM code^^ 
and typically employ around 600 processors simultane- 
ously 

Finally, we note that our CCM calculations based on 
the canted phase of the anisotropic kagome-latticc HAF 
shown in Fig. [TJb) do not assume for the s = \ model 
considered here that the canting angle <j> takes the same 
value 4>c\ = cos _1 (^-) as in the classical (s — > 00) 
case. Rather, calculations are first performed for an ar- 
bitrary choice of canting angle <j). We then minimize 
the corresponding LSUBm approximation for the energy 
EhS\JBm(4>) with respect to (f> to yield the corresponding 
approximation to the quantum canting angle ^lsub™- 
Generally (for m > 2) the minimization must be car- 
ried out computationally in an iterative procedure. Re- 
sults for the canting angle 0lsub«i are presented below 
in Sec. El 
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FIG. 2: (Color online) Ground-state energy per spin, E/Nk, 
of the spin-1/2 J1-J2 HAF on the anisotropic kagome lattice 
of Eq. ([TJ (with Ji = 1), using the LSUB4 approximation of 
the CCM with the ferrimagnetic canted model state shown 
in Fig. [lib), versus the canting angle <j> for various selected 
values of the anisotropy parameter Ji- For J 2 < 0.6 the min- 
imum is at <f> — (Neel' order) at this level of approximation, 
whereas for J2 > 0.6 the minimum occurs at 4> = 0lsub4 7^ 0, 
thereby providing evidence of a phase transition at J 2 ~ 0.6 
in this approximation. We show results for those values of 
4> for which the corresponding CCM equations having real 
solutions. 



IV. RESULTS 

We now present our CCM results for the anisotropic 
spin-1/2 J1-J2 HAF on the kagome lattice of Eq. ([TJ, 
where we use the canted ferrimagnetic state shown in 
Fig. IUb) as the CCM model (or reference) state |<I>) in 
the representations of Eqs. ([2]) and ([3]) of the exact GS kct 
and bra wave functions. Without loss of generality, but 
simply to set the energy scale, we henceforth set J\ = 1. 
(Equivalently, when more convenient to do so, we quote 
results in terms of the ratio n = J2/J1, assuming always 
that J\ > 0.) Wc first show in Fig. [2] the GS energy per 
spin, E/Nk-, as a function of the canting angle (j>. Al- 
though results are shown at the LSUB4 level of approx- 
imation results for other LSUBm levels are qualitatively 
similar. Curves such as those shown in Fig. [5] show that 
at this LSUB4 level of approximation, with J\ = 1, the 
GS energy is minimized at (f> = for J 2 < J 2 LSUB4 w 0.60 
and at a value <j> 7^ for J2 > J L B4 . Hence, these first 
results indicate the possibility of a slight shift of the crit- 
ical point to J2 = J2 1 = jLSUBoo k e t ween the quantum 
ferrimagnetic Neel' and canted phases, from the classi- 
cal value J2 = 0.5 when J\ = 1. The fact that Neel' 
order survives here, at least at finite orders of LSUBm 
approximation, into the regime where it would be clas- 
sically unstable against the formation of canted order is 
an example of a phenomenon that has been observed in 
many other magnetic systems, namely the tendency for 
quantum fluctuations themselves to favor collincar over 



FIG. 3: (Color online) The angle c^LSUBm that minimizes 
the GS energy i?LSUBm(0) of the spin-1/2 J1-J2 HAF on the 
anisotropic kagome lattice of Eq. (fTJ) (with J\ = 1), versus 
the anisotropy parameter J2. The LSUBm approximations 
with m = {2, 4, 6, 8}, using the ferrimagnetic canted state of 
Fig. 01b) as CCM model state, are shown. The corresponding 
classical result 4>c\ = cos^ 1 (^-) is also shown for comparison. 



noncollinear order. Nevertheless, we leave till later in 
this Section a discussion of cxtraplating these results for 
the critical point jLSUBm ^ Q ^ e m — > 00 limit, and hence 
also of a comparison of the quantum s = ^ case with its 
classical (,s — > 00) counterpart. 

In Fig. [3] we now show the canting angle 4>hSVBm that 
minimizes the GS energy i^LSUBm (</0 using the ferrimag- 
netic canted state of Fig. [lib) as CCM model state, at 
various CCM LSUBm levels, with m = {2,4,6,8}. Wc 
see from Fig. [3] that, at each LSUBm level shown, the 
canting angle </>lsubto that minimizes the corresponding 
estimate, i?LSUBm (</>), for the GS energy i?LSUBm(</>) of 
the model approaches zero smoothly, but with infinite 
slope, from the canted state (with 0LSUBm 7^ 0) side 
of the phase transition, as the anisotropy parameter J 2 
is reduced, to the corresponding estimate for the criti- 
cal point, jLSUB"^ and that it then remains zero for all 
■h < J^ sl!Bm on the Neel' state side of the transition. 
This behavior is completely analogous to that seen in 
the classical version of the model, also shown in Fig. [3] 
The evidence so far, therefore, is that both the classi- 
cal and spin-1/2 versions of the HAF on the anisotropic 
kagome lattice show second-order phase transitions be- 
tween ferrimagnetic states with collinear Neel' order and 
noncollinear canted order. Wc will sec below, however, 
that this picture changes when the magnetic order pa- 
rameter is considered too. 

We note first though that, by contrast, the correspond- 
ing behavior observed in the isotropic version (i.e., when 
J[ = Ji = 1) of the interpolating kagome-square HAF 
model (or J\-J' 1 -J2 model) of Fig. HJa), of which lat- 
ter model the present HAF model on the anisotropic 
kagome lattice is just the special case with J[ = 0, is 
quite different. Thus, in the spin-1/2 J\-J' 1 -J2 model 
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TABLE I: The critical value k£, at which the transition 
between the Neel' phase ((/> — 0) and the canted phase {<j> 7^ 0) 
occurs in the LSUBm approximation using the CCM with 
(Neel or) canted state as model state for the the spin-1/2 
J1-J2 HAF on the anisotropic kagome lattice of Eq. (0. 



Method 


,„LSUBm 


LSUB2 


0.685 


LSUB4 


0.600 


LSUB6 


0.568 


LSUB8 


0.559 
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with J[ = Ji = 1, it was observed 3 ^ that at each LSUBm 
level there is a finite jump in i^LSUBm at the correspond- 
ing LSUBm approximation for the phase transition at 
J 2 = jLSUBm between the Neel state (with L suBm = 0) 
and the canted state (with c/>LSUBm 7^ 0). This may be 
compared with the smooth behavior of the classical cant- 
ing angle </> c i = cos~ 1 (^) for this case for J 2 > Jf = 1 
for that model. For this latter Jx-J'^-Ji model in the 
isotropic case with J{ = J\ the evidence was that the 
phase transition between states with Neel and canted or- 
der was first-order for the spin-1/2 case compared with 
its second-order classical counterpart. On the other hand 
one should note that for that case^ we could not com- 
pletely rule out the possibility that as m — > 00, with in- 
creasing level of LSUBm approximation, the phase tran- 
sition at k — K C1 = K L SUBo ° becomes of second-order 
type, although a weakly first-order one seemed more 
likely on the basis of the available numerical evidence. 

Returning now to the question of estimating the phase 
transition point at k = k Ci in the present model, we note 
that previous empirical experienc o 36 ' 99 shows that the 
LSUBm estimates at k = K L s UBm fit well to an extrapo- 
lation scheme K^fUBm _ K LSUBoo +cm -i^ For p resen t 
anisotropic spin-1/2 J\—J<i model on the kagome lattice of 
Eq. , our phase transition estimates for k Ci = k lsuboo 
arc shown in Table fl] Using the above extrapolation 
scheme and the whole data set m = {2, 4, 6, 8} gives 
the estimate k Ci = 0.514 ± 0.003, while the correspond- 
ing estimate from using the data set m = {4, 6, 8} is 
n Cl = 0.515 ± 0.007. In both cases the quoted error esti- 
mates are simply the standard deviations from the asso- 
ciated least-squares fits. Our best estimate from combin- 
ing all of our results is k Ci — 0.515 ± 0.015, which may 
be compared with the corresponding classical value of 
K c i = 0.5. Clearly it is not excluded that k Ci — k c \, such 
that the transition point above which collinear Neel' or- 
der disappers occurs at exactly the same value n = 0.5 of 
the anisotropy parameter for both the extreme quantum 
and classical limiting cases of the spin quantum number 
s. We also discuss the nature of the phase transition at 
k Ci for the spin-1/2 model in more detail below. 

We note from Fig. [3] that as J 2 — > 00 the canting angle 
4> —> ^tt faster than does the classical analog <j> c \. Wc 
also note that for the special case J2 = 1 (or, equiva- 



FIG. 4: (Color online) Ground-state energy per spin, E/Nk, 
versus J2 for the spin-1/2 J1-J2 HAF on the anisotropic 
kagome lattice of Eq.{l]) (with Ji = 1), using the generic 
ferrimagnetic canted model state shown in Fig. [Tfb) as CCM 
model state, and with the canting angle <j> = 0LSUBm cho- 
sen to minimize the corresponding LSUBm estimate for the 
energy, -Elsubiti (<j>) , at each value of J2. The CCM LSUBm 
results with m = {2, 4, 6, 8} are shown, together with the cor- 
responding extrapolated LSUBoo result from Eq. 



lently, k = 1) of the isotropic kagome lattice, the CCM 
LSUBm estimates for the canting angle (j> take the value 
'/'LSUBm = 7^ for all values of m, as expected by symme- 
try, and exactly as in the classical version of the model. 

In Fig. [4] wc show our CCM results for the GS en- 
ergy per spin, E/Nk, as a function of J2, for the present 
anisotropic spin-1/2 J\- J2 model on the kagome lattice of 
Eq. ©. The CCM LSUBm results with m = {2,4,6,8} 
and the corresponding extrapolated LSUBoo results ob- 
tained from Eq. ((9|) are shown. As explained previously, 
in each LSUBm approximation we choose the canting an- 
gle <j) — </>LSUBm for each separate value of the parameter 
J2 that minimizes the corresponding CCM estimate for 
the energy, -E L suBm(»- 

At the isotropic kagome point (i.e., when J2 = Ji = 1) 
our present best estimate for the GS energy per spin, 
based on the extrapolation with the data set m = 
{4,6,8}, is E/Nk « -0.4352. This is lower than 
two recent rigorous upper bounds i 37 ' 38 Thus, Evenbly 
and Vidal^i used the multiscale entanglement renor- 
malization ansatz to evaluate exactly (up to floating 
point round-off errors) the energy of a wave function of 
the so-called honeycomb valence-bond crystalline type 
(with a 36-sitc unit cell), to give the rigorous bound 
E/Nk < —0.4322. Similarly, using a simple cluster 
product state for the infinite kagome lattice based on 
a fundamental cluster of 576 sites, for which the inte- 
rior of the cluster has the uniform valence-bond pat- 
terning expected of a spin-liquid state, Yan, Huse, and 
White 3 ^ have recently given an improved rigorous upper 
bound of E/N K < -0.4332. The same authors^ also 
use a large-scale density-matrix renormalization group 
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(DMRG) technique to provide what is certainly one of the 
most accurate estimates currently available for the en- 
ergy per site of the isotropic kagome-lattice HAF, namely 
E/N K = -0.4379 ± 0.0003. This estimate is itself con- 
sistent with the best available large-scale Lanczos exact 
diagonalization (ED) results for finite clusters of up to 
TV = 42 sites4£. It is also in excellent agreement with 
that from another very recent large-scale DMRG study^S 
namely E/N K = -0.4386 ± 0.0005. Our own present 
best estimate, cited above, is clearly below both of the 
rigorous upper bounds and in good agreement with the 
DMRG results. 

Our extrapolated result for the energy is also in very 
good agreement with previous CCM estimates^ for the 
spin- 1/2 HAF on the isotropic kagome lattice that used 
the kagome geometry itself to define the fundamental 
clusters of the LSUB?n configurations rather than the 
triangular lattice geometry of Fig. [lja) that we use here, 
as discussed above in Sec. IIIIl in which the Bi sites are 
retained. Thus, using the kagome-lattice geometry the 
extrapolated result E/Nk ~ —0.4357 was found with 
the LSUBm set m = {4,5,6,7,8,9,10} and E/N K « 
-0.4372 with the LSUBm set m = {6, 7, 8, 9, 10}4£ 

Figure [4] shows weak signals of the phase transition 
point (at k — Kci = ftLSUBoo) m Q f tne LSITB?™ 

curves, where a discontinuity in the first derivative of 
the energy is observed at the corresponding value k = 
K^ i SUBm . As usual the transition is seen more clearly in 
the behavior of the average local on-site magnetization, 
Mk = — -j^— Yli=i( s i)> where the sum is taken over all 



Nk sites of the kagome lattice and where again the spins 
are defined in the local, rotated spin axes in which all 
spins in the CCM model state point in the negative z- 
dircction. 

Thus, in Fig. [3] we show the magnetic order parame- 
ter, Mk, as a function of J 2 , for the present anisotropic 
spin- 1/2 J1—J2 model on the kagome lattice of Eq. (fTJ. 
CCM results are shown for LSUBm approximations with 
m = {2,4,6,8}, together with various LSUBoo extrapo- 
lations. As discussed in Sec. IIII| for the strongly frus- 
trated regime in which the canted state is the stable 
ground state (i.e., for J 2 > J2 1 ) we use the well-tested 
and established scheme of Eq. (|10p , whereas for the less 
frustrated regime in which the Neel' state is the stable 
ground state (i.e., for J 2 < J| L ) we use the scheme of 
Eq. pT|) . Since for the Neel' state the LSUBm results 
converge (with increasing values of m) much faster than 
those for the canted state, as can clearly be seen from 
Fig. El it is evident that the use of these different schemes 
for the two regimes is justified. Since the approximate 
transition point at J2 = jL SUBm (when J\ = 1) between 
the two phases depends slightly on the CCM truncation 
index m, as has already been noted above, and as can be 
seen clearly in Fig. [5] (and, more explicitly in Table HJ), it 
is clear that the region very near the transition is inher- 
ently difficult to extrapolate accurately. 

To illustrate the sensitivity of our extrapolations to the 
approximations used we show in Fig. [5] the corrcspond- 
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FIG. 5: (Color online) Ground-state magnetic order pa- 
rameter, Mk, versus J2 for the spin-1/2 J1-J2 HAF on the 
anisotropic kagome lattice of Eq. fT]) (with J\ = 1), using the 
generic ferrimagnetic canted model state shown in Fig. QJb) 
as CCM model state, and with the canting angle (f> — ^suBm 
chosen to minimize the corresponding LSUBm estimate for 
the energy, -BLSUBm(0), at each value of J2. The CCM 
LSUB771 results with m = {2, 4, 6, 8} are shown, together with 
the corresponding extrapolated LSUBoo results using both 
data sets in = {2, 4, 6, 8} and m = {2, 4, 6} for comparison 
purposes. Extrapolated results are calculated using Eq. (|11[) 
for the Neel' state and Eq. (|10|l for the canted state, respec- 
tively; and the (grey) dashed vertical line at J2 = 0.515 repre- 
sents our best estimate for the termination point above which 
collinear Neel' order disappears, at J2 = J^ 1 = j^suboo^ ag 
discussed in the text. 



ing extrapolations in the two regimes using both the data 
sets m = {2,4,6,8} and m = {2,4,6}. Since the (most 
accurate) LSUB8 scheme is computationally expensive, 
results are shown only at the limited set of J2 values in- 
dicated by the symbols in Fig. [5] It is very encouraging 
that the extrapolated curves for Mk are very steep near 
J2 1 and that they become steeper still as higher LSUBm 
approximations are included. For example, the extrap- 
olated LSUBoo curve for Mk obtained from the set 
m = {2,4,6} becomes zero at J 2 rj 0.67 (with J\ = 1), 
while that obtained from the set m = {2, 4, 6, 8} becomes 
zero at J2 ps 0.63 (with J\ = 1). It seems reasonable to 
assume that the actual LSUBoo curve will become verti- 
cal at the point J 2 = J C 2 X = J^ SVBoc ^ as is seen in Fig. [5] 
by the proximity of the extrapolations to the vertical line 
at J2 = 0.515, which represents our best estimate for the 
transition point between the Neel' and canted states. 

Figure [S] shows the existence of a clear window in the 
parameter k in which Mk < 0, and hence in which 
the canted order present in the model state has van- 
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ished. This (paramagnetic) region includes the point 
K = 1 corresponding to the isotropic kagome HAF, and it 
seems reasonable to assume that the phase present in this 
regime is the same as the paramagnetic GS phase of the 
isotropic model. As discussed above, the lower boundary 
of this window seems to coincide with the point n = n Cl 
above which the collinear Neel' order disappears. 

Thus, the evidence from the magnetization data shown 
in Fig. [5] now indicates that the transition at k = n ci 
is actually between the Neel' and paramagnetic phases, 
rather than between the Neel' and canted phases as in 
the classical case at the corresponding value k = n c \. 
Nevertheless, we cannot completely exclude the possibil- 
ity of a very narrow strip of canted phase between the 
Neel' and paramagnetic phases confined to the region 
0.5 < k < 0.6. 

We denote by n C2 the corresponding critical value of k 
that marks the upper boundary of the window in which 
Mk < 0, and that hence marks the transition between 
the paramagnetic and canted phases. From Fig. [5] we 
find estimates k C2 ps 1.83 ± 0.02 from the extrapolated 
LSUBoo curve using the data set m = {2,4,6,8} and 
k C2 ps 1.80 ±0.02 from the corresponding curve using the 
data set m = {2, 4, 6}. 

In Fig. [6] we also show, in various approximations, the 
separate average on-site magnetizations, Ma and Mb 2 , 
on the A sites and the B2 sites respectively of Fig.[IJb), as 
well as their average value on the kagome lattice, Mk = 
§ Ma + \M-q 2 . We recall that the A sites are connected 
by two Ji bonds and two J 2 bonds, while the B 2 sites 
are connected by four J\ bonds. In particular, we note 
that as J 2 — > 00 (with J\ = 1), the model reduces to one 
of independent linear HAF chains of alternating Ai and 
A2 sites. In this limit our extrapolated CCM result for 
the energy using the ferrimagnetic canted state as model 
state (and see Fig.gjl is E/N K ps -0.2954 J 2 . Since in this 
limit the B 2 sites become irrelevant, we may re-express 
the result in terms of the number, N c = fA^-, of A sites 
that form the independent ID chains in this limit, as 
E/N c ps — 0.4431</ 2 . This compares extremely well with 
the exact result for the ID HAF, obtained from the Bethe 
ansatz solution, E/(N C J 2 ) = \ - In 2 ps -0.443147. 

As J 2 -> 00 (with Ji = 1) all of our CCM LSUBm ap- 
proximations give Mb 2 —> 0.5, as expected. The LSUB4 
result for Mb 2 shown in Fig. [5] is typical of the behav- 
ior of the entire LSUBm set. By contrast, in the same 
limit, the individual LSUBjti approximations for Ma ap- 
proach different constant (positive) values. The extrap- 
olated LSUBoo result for Ma, shown in Fig. [6l however, 
approaches a value very close to zero, again fully consis- 
tent with the exact behavior of ID HAF chains. 

Our results so far have indicated the presence of a para- 
magnetic state in the range k Ci < k < n C2 between the 
quasiclassical states with Neel' order (for k < k Ci ) and 
canted order (for k > n C2 ). We have seen too that as 
k — > 00 the model reduces consistently to the limit of un- 
coupled isotropic HAF ID chains. There remains still the 
question raised in Sec. [H] as to whether the canted state 
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FIG. 6: (Color online) Various ground-state magnetic or- 
der parameters versus J2 for the spin-1/2 J\—J% HAF on 
the anisotropic kagome lattice of Eq. (fl]) (with Ji = 1), us- 
ing the generic ferrimagnetic canted model state shown in 
Fig. [Ub) as CCM model state, and with the canting angle 
4> = <^LSUBm chosen to minimize the corresponding LSUBm 
estimate for the energy, E^svBm(4>), at each value of J2. Re- 
sults are shown for the average local on-site magnetizations, 
Ma and Mb 2 , on the A sites and the B2 sites respectively 
of Fig. [Ub), as well as for their average value on the kagome 
lattice, M K = §M A + |Mb 2 - LSUB4 results are shown for 
Mb 2 , and extrapolated LSUBoo results are shown for Mk us- 
ing both data sets m = {2, 4, 6, 8} and m — {2, 4, 6}, and for 
Ma using the data set m = {2, 4, 6}. Extrapolated results are 
calculated using Eq. for the Neel' state and Eq. (|10|l for 
the canted state, respectively; and the (grey) dashed vertical 
line at J2 = 0.515 represents our best estimate for the termi- 
nation point above which collinear Neel' order disappears, at 
J 2 = J2 1 = jL SUBoo 5 as discussed in the text. 



illustrated in Fig. [Ub) remains the stable GS phase all 
the way out to k — > 00, where the canting angle <p — > 90°, 
or whether there might exist a further phase transition at 
a value k — k C3 > k C2 to some other phase. One such pos- 
sibility is the semi-striped state illustrated in Fig. [IJc) . 
We have argued that such a phase might be stabilized due 
to the quantum fluctuations possibly lifting the infinite 
degeneracy, which exists in the classical counterpart be- 
tween the orientation of the antiferromagnetically aligned 
spins on the A sites and the orientation of the ferromag- 
netically aligned spins on the B 2 sites (that become de- 
coupled from those on the A sites in this limit), by the 
order-by-disordcr mechanism. We consider this possibil- 
ity further in Sec. El where we also discuss our results 
and compare them with those of others. 



V. DISCUSSION AND CONCLUSIONS 

We have used the CCM to investigate the effects of 
quantum fluctuations on the zero-temperature GS prop- 
erties and phase diagram of the spin-1/2 Ji-J 2 HAF 
on the anisotropic kagome lattice of Eq. (TTJ), and as il- 
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lustratcd in Fig. [TJb). The system contains spatially 
anisotropic NN exchange couplings on the kagome net, 
with coupling J 2 = kJ% > in one of the three equiva- 
lent spatial directions of the lattice and coupling J\ = 1 
along the other two directions. The model has only two 
classical GS phases. For k < k c \ = i the GS phase has 
collinear fcrrimagnetic Neel' order, in which the J2-chain 
spins on the A sites of the lattice are aligned in one di- 
rection and the middle spins on the remaining B2 sites 
are aligned in the opposite direction. At k = k c \ = i the 
classical GS configuration changes from the essentially 
unique collinear ferrimagnetic Neel' state to an infinite 
ensemble of degenerate states. These include a still infi- 
nite number of coplanar canted ferrimagnetic states that 
are expected to be selected from among the rest by ther- 
mal or quantum fluctuations to comprise the stable GS 
phase for all k > n c \ = i. 

For a given value of n > k c \ all of these states are char- 
acterized by a canting angle such that on each A1A2B2 
triangular plaqucttc of the kagome net the two J2-chain 
spins (i.e., those on the Ai and A 2 sites) form angles 
(ir ± </>) with respect to the middle spin (i.e., that on 
the B2 site). Clearly, the Neel' state is just the special 
case with <fi = 0. Each triangular plaqucttc in the en- 
semble of coplanar states also carries a chirality variable, 
X = ±1, defined to be the direction (anticlockwise or 
clockwise, respectively) in which the spins rotate as one 
traverses the plaquette in the positive (anticlockwise) di- 
rection. The different degenerate coplanar canted states 
then correspond to different ways of assigning chiralitics 
to the individual plaquettes. 

In this paper we have used for the spin-1/2 model 
the generic ferrimagnetic canted model state shown in 
Fig. [Ub) as our CCM model state, which is that mem- 
ber of the classically degenerate ensemble in which all of 
the triangular plaquettes have the same (here positive) 
chirality. This state thus corresponds to the q = state 
of the isotropic (i.e., when J2 = 1) kagome-lattice HAF, 
in which </> = ^ by symmetry. At each LSUBm level 
of approximation we have chosen the value of the cant- 
ing angle (j) that minimizes the corresponding LSUBm 
estimate for the GS energy. We note again parenthet- 
ically that in a previous recent CCM analysis 4 ^ of the 
isotropic kagome HAF in which both the q = state and 
the V3 x \f?> state (that corresponds to that member of 
the classically degenerate ensemble in which the chirali- 
ties alternate, such that triangular plaquettes joined by 
a vertex have opposite values of x) > it was found that for 
the extreme quantum case considered here, with s = i, 
the q — state is energetically favored over the y/3 x \/3 
state, while for any s > i the \J~S x \/3 state is selected 
over the q = state. 

We found here that the canting angle 0LSUBm that 
minimizes the GS energy i?LSUBm(0) for the spin-1/2 
model, at a given LSUBm level, becomes nonzero for 
values of the anisotropy parameter k > K^ i SUBm , where 
K ci SUBm i s generally somewhat higher than the corre- 
sponding classical value n c \ = ^ , as shown in Tabic |TJ On 



the other hand these critical values converge to an ex- 
trapolated value k Ci = k l suBo ° that is close to « c i = -|. 
Indeed our best estimate is k Ci = 0.515 ± 0.015. Al- 
though they cannot entirely exclude the possibility of a 
small regime of canted phase in a very narrow strip im- 
mediately above k Ci , our corresponding CCM LSUBm 
results for the magnetization shown in Figs. [5] and |6] give 
compelling evidence, however, that the transition for the 
spin-1/2 model at k = k Ci from a GS phase with Neel' 
order is not to one with canted order, as in the corre- 
sponding classical model at K = n c \ = 0.5, but rather 
to a paramagnetic state with no canted order. Our ev- 
idence is that this paramagnetic GS phase persists over 
the anisotropy range k Ci < k < n C2 , before the canted 
state becomes the stable GS phase for k > k C2 . Our best 
estimate for this upper critical point of the paramagnetic 
phase is k C2 = 1.82 ± 0.03. 

Since the isotropic kagomc-lattice point, k = 1, is con- 
tained within the parameter range k Ci < k < k C2 of this 
paramagnetic phase, the natural conclusion is that this 
phase of the anisotropic model shares the same order as 
the GS phase of the isotropic model. As we discussed in 
Sec. [H the isotropic spin-1/2 HAF on the kagome lattice 
has been greatly studied in the past. The most direct 
results, namely those from the exact diagonalization of 
finite lattices ; 14 ' 16 seem to provide strong evidence for a 
spin-liquid GS phase. This conclusion is supported by the 
results of block-spin approache d 15 ' 18 and by those from 
various other studies toOf 5 -^^^^^^ 

On the other hand, conflicting results have been found 
by other author a 3 ' 21 ! 22 ' 24 ' 27 ' 30 ' 37 who have proposed var- 
ious valence-bond solid states as the GS phase of the 
isotropic HAF on the kagome lattice. A detailed compar- 
ison of the exact spectrum of a 36-sitc finite lattice sample 
of the isotropic kagome HAF against the excitation spec- 
tra allowed by the symmetries of the various proposed 
valence-bond crystal states has, however, cast very strong 
doubts on their validity as stable GS phases.— Over the 
past year or so this muddled and confused picture of the 
nature of the GS phase of the isotropic spin-1/2 HAF on 
the kagome lattice has been seemingly resolved in favor 
of a topological spin liquid. 

In particular, as indicated in Sec. U the results of two 
independent and very recent large-scale DMRG studies 
of the spin-1/2 isotropic HAF on the kagome lattic o 49 ' 50 
have provided compelling positive evidence that this GS 
phase is a topological quantum spin liquid. The simplest 
such state that preserves all symmetries is the Z2 spin 
liquid, and by explicitly calculating the topological en- 
tanglement entropy both recent DMRG studies provide 
strong positive evidence that the spin liquid state does 
indeed have Z2 topological order. 

It has not been our aim here to investigate the order 
properties of the paramagnetic GS phase of antifcrro- 
magnetically coupled 5 = 5 spins on the infinite kagome 
lattice, but rather to investigate the stability of the phase 
as the lattice is spatially distorted. Clearly, however, it is 
natural to expect that over the entire range k Ci < k < k C2 
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in which the paramagnetic phase persists (as manifested 
here by a negative, and hence unphysical, value of the 
calculated local magnetic order parameter) it remains a 
spin liquid with the same topological order. 

The spin- 1/2 HAF on the spatially anisotropic kagome 
lattice has also been studied by several other authors 
recently using a variety of techniques. These have 
included large-TV expansions of the Sp(iV)-symmetric 
generalization of the actual SU(2) model , 77 ' 78 a block- 
spin perturbation approach to the trimerized kagome 
lattice^ 8 - various semiclassical calculations (appropriate 
to the limit of large spin quantum number s) that in- 
clude (a) studying an effective chirality Hamiltonian de- 
rived from a low-temperature classical nonlinear spin- 
wave expansion^ and (b) keeping terms of order 1/n 
in the large-n limit of the O(n) generalization of the 
classical 0(3) model together with a high-temperature 
expansion,— field-theoretical techniques appropriate to 
quantum critical systems in one dimension (and which 
are hence appropriate here for the case re = J2/ Ji 3> 1 of 
weakly coupled chains)^ and a renormalization-group 
analysis in the same quasi- ID limit but now also in the 
presence of a Dzyaloshinskii-Moriya interaction^ 

Since many of these calculations employ perturbation 
theories of one kind or another in some "artificial" small 
parameter, direct comparison is difficult. Thus, for ex- 
ample, in the largc-iV Sp(7V) expansion, the effective 
smallncss parameter is a = n^/N, where rib is the num- 
ber of bosons on each site. While in the physical SU(2) 
model (which corresponds to the case N ~ 1) we have 
a = 2s, the comparison in the large-iV limit actually 
studied is lost. Yavors'kii et a/i 7 - 8 - argue that the value 
of a that corresponds to the s = \ under study must 
be somewhat less than 0.5. Similarly, in the block-spin 
perturbation approach of Yavors'kii et alS^ the assump- 
tion is made that the kagome lattice is trimerized such 
that the spins on the downward-pointing triangles, say, 
are strongly coupled whereas the couplings on the bonds 
of the upward-pointing triangles are weaker by a factor 
7. The approximate GS phases of this trimerized model 
are then studied in different regimes of the anisotropy pa- 
rameter k in a perturbation expansion with respect to 7, 
while the physical model corresponds to the case 7 = 1. 

Complementary to such essentially perturbative stud- 
ies have been various more controlled analyses of the 
quasi-lD limit of the mode h 80 ' 82 The latter studies, by 
their nature of focussing on the large-anisotropy (k » 1) 
limit, also lose sight of the intermediate paramagnetic 
(spin- liquid) phase that has been the focus of the present 
study. For example, Zyuzin et alJ& explicitly state that 
they do not find a spin-liquid ground state in any regime 
that they study. 

As we have seen, much of the previous work on the 
spin- 1/2 HAF on the spatially anisotropic kagome lat- 
tice has approached the quantum limit only very indi- 
rectly, either from the classical side or in such slave par- 
ticle approaches as the Schwinger boson technique ap- 
plied in the large-A^ Sp(N) approach. The only direct 



s = -| approach seems to be a small-scale ED study of 
up to Nk = 24 spins (in a 4 x 2 unit cell arrangement) 
with periodic boundary conditions.— The numerical ev- 
idence from the ED study seems to indicate very clearly 
that for values of the anisotropy parameter ft < 0.5 the 
GS phase has nonzero total spin. Indeed, its value for 
the few finite-size lattices studied is precisely what is ex- 
pected for the classical collinear Neel' fcrrimagnctic state, 
namely S'tot = ^Nks, thereby agreeing fully with our 
own findings. By contrast, the numerical evidence for 
values k > 0.5 is far less clear. While the evidence seems 
to be that for all lattice sizes up to 24 sites the GS phase 
is a spin singlet for all values k > 0.5, there is a clear ten- 
dency for a state of nonzero spin at the T point to drop 
in energy on moving away from the isotropic point re = 1, 
perhaps indicating a tendency to develop a net moment 
again. Nevertheless, the evidence from such small-scale 
ED studies seems to leave completely open the nature of 
the GS phase for k > 0.5. 

The two semiclassical approaches of Wang et al£2- also 
concur that for k < 0.5 the GS phase is the collinear 
Neel' fcrrimagnetic state. For the case re > 1 both ap- 
proaches also indicate a canted fcrrimagnetic state of the 
classical type, but where the infinitely-degenerate mani- 
fold of coplanar states is lifted by the order-by-disorder 
mechanism, which now seems to favor the so-called chi- 
rality stripe state as the GS phase, in which all spins on 
the interstitial B2 sites are ferromagnctically aligned, and 
where pairs of triangles on the kagome lattice that share 
either an Ai or A2 vertex have the same value of the chi- 
rality parameter \i while pairs sharing a B2 vertex have 
opposite values of \. (We note parenthetically that this 
state is the only other state, apart from the q = state 
of Fig. [ljb), that is allowed by the chirality constraints 
in the general case k / 1 and in which all of the B2- 
spins are aligned parallel to one another.) By contrast, 
in the intermediate regime 0.5 < re < 1, the semiclassical 
models do not give clear indications of which GS order- 
ing is favored. Thus, while the y/2> x a/JS order seems to 
be favored in the spin-wave expansion for the isotropic 
case re = 1, the comparable analysis for the k < 1 case 
seems to depend very sensitively on both the choice of 
unphysical parameters and on the number of chirality- 
chirality couplings included in the analysis. By contrast, 
the large-n saddle-point solution of the 0(n) generaliza- 
tion of the classical 0(3) model seems to favor the q = 
GS ordering of spins for the case 0.5 < re < 1. 

The l&rge-N Sp(7V) expansion analysis of the 
mode l 77 ' 78 indicates that the actual spin-1/2 HAF on the 
spatially anisotropic kagome lattice has a GS phase with 
collinear Neel' fcrrimagnetic order for small values of the 
anisotropy parameter k, which gives way at larger val- 
ues of k to an incommensurate (spiral) GS phase with no 
LRO, in general agreement with our findings. In turn, 
as k is increased further this GS phase then gives way 
to a phase in which the chains arc completely decoupled, 
while the interstitial spins (i.e., those on B2 sites) show 
some short-range spin-spin correlations. For reasons al- 
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ready noted above this Sp(7V) analysis is unable to give 
quantitative estimates for the corresponding two critical 
values of k for the actual spin-1/2 SU(2) model. 

Finally, we note that the canted ferrimagnetic phase 
that we have found to be the stable GS phase for k > k C2 , 
after the disappearance of the paramagnetic (spin-liquid) 
phase is not likely to remain the GS phase for sufficiently 
large values of k, as we have already noted in Sec. [TT1 
since in the limit n — > oo the canting angle (/>—>■ ^, and 
the chain spins become perpendicular to the interstitial 
spins. Since in this limit the relative orientation of the 
chain spins and the interstitial spins becomes irrelevant, 
and since quantum fluctuations generally prefer collinear 
spin configurations in such a degenerate situation, we 
thus expect a third phase transition at a value k = n C3 
to a phase that eventually becomes the decoupled ID 
HAF chain phase in the asymptotic limit k — > oo. The 
precise nature of this fourth phase is by no means settled. 

We have alluded in Sec. [TT] to one such candidate be- 
ing the collinear ferrimagnetic semi-striped phase shown 
in Fig.QTc). We are currently investigating whether this 
phase might become energetically favored at sufficiently 
large values of k. We intend to report on this in a sepa- 
rate future paper. 

An alternative candidate state for the large-K phase 
has been suggested by Yavors'kii et al~ from their block- 
spin trimerized version of the model. In the subsequent 
small-7 perturbative limit they find a tentative ground 
state in the large-anisotropy case (k ^> 1) that is a 
collinearly ordered antiferromagnet in which the inter- 
stitial (B2) spins are Neel-ordered and the spins on the 
A-chains form singlet dimcrs (i.e., the spins on each Ai- 
A2 NN pair of sites on, say, each downward-pointing tri- 
angle on the kagome lattice form a spin-singlet state). 

Yet another analysis of the large-K limit, by Schnyder 
et a/.,— suggests that all of the spins order with a (gener- 
ally noncoplanar) configuration in which the interstitial 



spins on B2 sites and the chain spins on A-sites each sep- 
arately form predominantly coplanar spirals with a wave 
vector (g,0), but with a reduced [0(1/ n)} static moment 
on the J2-coupled chains. These authors find that the 
chain spins are weakly canted out of the plane, with the 
[0(1/ k 2 )] normal components being ordered in an an- 
tiferromagnctic fashion. While their analysis could not 
determine q reliably, it is expected that g« 1 and, in- 
deed, the authors suggest that q = is a real possibility, 
in which case the state becomes coplanar. Nevertheless, 
even this state differs from both the chirality stripe state 
considered by Wang et al3^ and the dimerized state con- 
sidered by Yavors'kii et alZ^, although there are some 
similarities with each. 

It is clear that the large-K limit of the spin-1/2 HAF 
on the anisotropic kagome lattice is far from settled. Our 
own work presented here has mainly been concerned to 
investigate the stability with respect to anisotropy k of 
the spin-liquid state that has convincingly been found in 
recent work to be the stable GS phase of the isotropic 
(k = 1) model. Nevertheless, we hope to return in the 
future to the quite separate question of whether or not 
there is a further transition at some value k — n C3 > n C2 
of the anisotropy parameter, from the canted state dis- 
cussed here to some other state with or without collinear 
order. 
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